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Abstract. The blue-shifted broad emission lines and/or broad absorption lines seen 
in many luminous quasars are striking evidence for a broad line region in which radia- 
tion driving plays an important role. We consider the case for a similar role for radiation 
driving beyond the dust sublimation radius by focussing on the infrared regime where 
the relationship between luminosity and the prominence of the 3-5 /mi bump may be 
*-£h ■ key. T o investigate this further, we apply the 3D hydrodynamic wind model of lEverettl 

*~"! (l2005h to predict the infrared spectral energy distributions of quasars. The presence 

of the 3-5 yum bump and strong, broad silicate features can be reproduced with this 
dynamical wind model when radiation driving on dust is taken into account. 



1. Introduction 



The torus was introduced by lAntonucci & Milled dl985 ) to explain the revelation from 



spectropolarimetry that Type 2 (narrow-line) Seyfert galaxies can look like Type 1 
(broad-line) Seyferts along different lines of sight. Since that discovery, the torus' 
primary purpose has been to obscure the broad-line region, but it also serves to ex- 
plain the near- and mid-infrared (IR) emission observed in the vast majority of quasars 
whose spectral energy distributions (SEDs) indicate dust heated to near sublimation 
temperatures, presumably by the accretion-disk continuum. 

Typically , the torus has been modeled as a static, axisymmetric structure (e.g., 
Pier &KrohM 1 19921) hat may have condensations (clumps or clouds) within it (e.g., 
Nenkova et al.ll2008l). A notable exception to this paradigm is the dusty wind model of 
Konigl & Kartje ( 1994), who started with a magneto-hydrodynamic (MHD) wind, and 



added in radiation driving on grains beyond the dust sublimation radius. They were 
motivated in part by the large scale-height of the obscurer required by the observed 
Seyfert 1 to Seyfert 2 ratio, the inevitability of radiation driving in the presence of such a 
strong sourc e of UV photons , and the apparent lumino sity-dependent covering fr action 
of the torus dLawrencelll99ll) . Since then others (e.g. jElitzur & Shlosman 2006) have 



also incorporated some wind component into their m odels, and/or the effects of IR 
radiation-pressure in inflating the dusty medium (e.g.. iDorodnitsyn et al.ll20lT ). The 



questions remains: how to test for the existence of this dusty wind? 

In the UV, wind signatures are clearly evident in the broad resonance lines of high 
ionization species such as C iv. W hile most apparent in classic P-Cygni-like broad 
absorption lines (e.g., iLy ndsl [l967h . the asymmetries in the same broad emission lines 
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Figure 1 . A composite IR spectrum made from ~ 50 luminous SDSS quasars from 
the Spitzer IRS archive. Prominent emission features are labeled; the horizontal bars 
indicate the approximate location of the broad silicate emission bumps. From Hill et 
al. (in preparation). 



are also evidence for winds (e.g.. lLeighlyll2004l ; lKruczek et al.ll201ll) . The most natural 
mechanism for accounting for the high velocities (~ 10 4 km s~ ! ) of the outflows is 
resonance-line driving by UV continuum photons. 

Unlike the clear wind signatures in the UV, quasar IR spectra are characterized 
by broad features, such as the silicate emission bumps from the vibrational stretching 
and bending modes of the grains at 10 and 18 /mi, respectively (see Figured]). The 
widths and peaks of the silicate bumps are affected by the grain size distribution, and 
so extracting information on the velocity of the emitting grains is impossible given the 
uncertainties in the inherent shape of the bumps. Narrow atomic features in the IR are 
from forbidden lines such as [O iv] and [Ne v] generated in low density, photoionized 
gas in the inner parts of the narrow-line region. Clear spectral signatures of a dusty 
wind are therefore elusive. 



2. Luminosity Dependence of the 3-5 jum Bump 



In the UV, the maximum possible ter minal velocity of outflows seen in absorpti on ap- 
pears to be set by the UV luminosity dLaor & Brandtll2002l ; rGanguly et al.ll2007l) . This 
makes sense for a radiation-driven UV outflow; a high UV photon density has the ability 
to transfer more momentum to the ionized wind. While this particular mechanism is not 
relevant in the IR, spectral features in this regime whose strength depends on luminosity 
are potentially fruitful for looking for evidence of the wind from the dust. One such fea- 
ture is the 3-5 /urn bump, which can be successfully modeled as a b lackbody with a tem- 
perature of 1200 — 1500 K (e.g.. lBarvainis|[T987l ; lDeo et al.ll201lh . These temperatures 
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are interesting because they are near the grain sublimation te mperature of graphi t es, bu t 
higher than the grain sublimation temperature of silicates. lEdelson & Mallear] (1 19861) 
first n oted that this feature becomes stronger at higher luminosities. |Gallagher et al. 
(120071) confirmed this trend with a much larger sample of SDSS quasars with Spitzer 
IRAC + MIPs photometry by searching for convex spectra between 1 and 8 jira; the 
amount of spectral curvature was significantly correlated with luminosity. 

The increasing prominence of the 3-5 //m bump can be qualitatively understood 
in a wind paradigm, because more luminous quasars will have more radial IR outflows. 
Assuming an accretion disk is the source of the optical through X-ray continuum, the 
dust-driving continuum will always be interior to (and significantly more compact than) 
the dusty wind, and thus act in the radial direction. In this scenario, the hottest dust is 
visible from a larger range of inclination angles in luminous objects. In the UV, an 
ionized wind launched from the accretion disk at small radii (~ 10 16 cm) could have a 
significant vertical component to its acceleration. This geometric effect can naturally 
account for the larger observed fraction of BAL quasars (seen when looking through 
the UV wind) compared to type 2 quasars (seen when the broad-line region is blocked 
by thick, dusty material). 



3. The MHD Dusty Wind Model 



To investigate the interplay of radiation driving and the IR SED, we have developed a 
dynamical model of the torus as a dusty wind launched by MHD forces and by radiation 
pressure from the accretion disk continuum. The dusty wind generated in this manner 
can cover a large fraction of the sky as seen from the central black hole. This model has 
the benefit of being more self-consistent than static torus models and including the im- 
portant physics of motion around the black hole and radiative acceleration. Ultimately, 
we wish to use this model to understand how the physical properties of dusty winds in 
quasars correlate with their observable spectral signatures in the IR. 

The mod el and its correspondi ng code has been expanded from its original form 
(described in lKonigl & Kartell 19941) bv lEverettl(l2005l) . where a comprehensive account 
of the model's components and key equations are described. We advance on lEverettl 
( 20051) by a dding the continuum opacity of ISM du st grains, as specified by the ISM 
dust model dMathis et al.ll 19771 : Ivan Hoof et al.ll2004l) in Cloudy dFerland et al.lll998h . 

The radiation pressure on the dust from the central source is very strong - even 
for cases where L/Ledd = 0.1, it is appr oximately 1 times the force required to un- 
bind dust from the gravitational potential (lEverett et al.ll2009h . Dust grains absorb radi- 
ally streaming photons originating from the accretion disk, and then cool by radiating 
isotropically. The force, felt by the dust particles due to conservation of momentum, 
feeds back on the wind structure by bending the magnetic wind from a vertical to more 
radial structure. The radiative force works in conjunction with the magneto-centrifugal 
forces to accelerate the wind flow and therefore modifies the structure of the outflow. 

After iterating the MHD+radiation pressure calculations to se t up the stru cture of 
the wind, we use the Monte Carlo radiative-transfer code MC3D dWol J 120031) to send 
rays through the wind and predict the observed IR SEDs. The input parameters for 
the wind such as Eddington ratio (L/Lndd), black hole mass (Mrh), and illu minating 
continuum (the empirical optical-to-X-ray composite of iRichards et al.i r2006. hereafter 
R06) are varied to investigate their effects on the output SED. 
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Figure 2. Left: IR to X-ray quasar SEDs as follows: the empirical R06 composite 
of optically luminous quasars (dotted); the input incident accretion disk continuum 
(dashed); the output IR SED for the characteristic model (red solid); the summed 
output IR SED and input incident accretion disk continuum (solid black). A char- 
acteristic model with a smaller dust-sublimation radius shows an SED (blue solid) 
peaking in the near-IR. As the dust-sublimation radius is increased the SED power 
decreases in the near-IR and the peak emission shifts to longer wavelengths; small 
dust sublimation radii are required to generate strong 3-5 //m bump emission. Right: 
The output IR SEDs generated by the wind models; different line styles are used to 
distinguish the observed inclination angle, i, with respect to the normal to the ac- 
cretion disk. The power and shape of the SED is not particularly sensitive to the 
angle, indicating that the wind is optically thin to IR radiation for A > 8 fim. The 
short wavelength emission increases with smaller inclination angle. This effect could 
contribute to relatively stronger 3-5 //m emission in more luminous objects because 
their more radial dusty winds reveal the hottest dust in a larger range of inclination 
angles. Figures from Keating et al. (2012). 
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4. Principal Results 

With input parameters appropriate for luminous quasars (luminosity of the incident 
continuum: Lt, i ~ 10 46 erg s" 1 ; L/L^dd = 0.1; Mbh = 10 8 M Q ; column density at the 
base of the wind: A^h,o = 10 25 cm -2 ; and observer's inclination angle from the disk 
normal: i = 60°), we can produce reasonable IR SEDs with approximately the right 
shape and luminosity (L « 10 43 - 10 45 erg s -1 ) as expected from the R06 composite 
of optically luminous SDSS quasars (see Figure [2]). This is a promising result for a 
relatively simple model given that we have not attempted any fitting. 

A benefit to our model is that we are able to see directly the effects of various 
physical parameters on the final IR SED, which ultimately will allow us to understand 
the physical properties of the torus itself. By determining which physical parameters 
have an observable effect on the IR SEDs, and narrowing them down so that we gen- 
erally reproduce the power expected, we have established a reasonable starting point 
from which we can expand and further refine our model. 



5. Future Work 

In the near future, we plan to implement a more sophisticated treatment of dust grains 
that takes into account the different sublimation temperatures of graphite and silicate. 
The luminosity dependence of emission from the hottest grains — graphites — is 
promising as a means of demonstrating the connection between radiative driving of 
grains, wind geometry, and the observed SED. We expect that radiation pressure must 
be important, because otherwise the luminosity dependence of the torus covering frac- 
tion and of the 3-5 /mi emission is hard to understand. 

Our ultimate goal is to generate a library of realistic, IR SEDs that can be com- 
pared with Spitzer IRS observations of quasars. We aim to determine (1) if reasonable 
input properties (e.g., the input continuum shape, L/L^dd, and Mbh) result in SEDs that 
match in detail those observed, and more ambitiously, (2) how to use the empirical IR 
SEDs to independently constrain the unknown physical properties of quasars. 
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